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Abstract 

We study the formation of domains in a continuous phase transition with a 
finite-temperature quench. The model treated is the <3? 4 theory in two spatial 
dimensions with global 0(2) symmetry. We investigate this using real-time 
thermal field theory, following Boyanovsky and collaborators, and find that 
domain sizes appear to be smaller than those produced in an instantaneous 
quench in the tree-level approximation. We also propose that a more physical 
picture emerges by examining the two-point functions which do not involve 
any cutoff on the short wavelength Goldstone modes. 



I. INTRODUCTION 

The formation of domains in phase transitions, and the subsequent (growth) dynamics 
of these domains is an important topic arising in a variety of contexts ranging from early 
universe cosmology to condensed- matter physics. As a result, domain formation has been 
much investigated from both the classical |l| and quantum mechanical |2| points of view. 

There are two principal mechanisms of domain formation, namely a) Bubble-Nucleation 
and b) Spinodal Decomposition. Strongly first order phase transitions typically proceed 
via bubbie-nucleation while spinodal decomposition is often observed in second order phase 
transitions. For a weakly first order phase transition, however, the idea that the phase 
transition is driven by the formation of bubbles is not so clear. In fact, this type of phase 
transition can be driven by spinodal decomposition as well. Nevertheless, domain-type 
structures appearin both mechanisms. 

The growth process of such domains at late times (also known as phase ordering in the 
literature) is relatively well-understood [0 while the study of the formation process of the 
defects associated with these domains has begun only recently 0. The initial size of the 
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domains formed in continuous phase transitions is of physical interest. For phase transitions 
that generate topological defects this domain size also determines the initial defect density. 

The earliest attempt to determine the initial domain size, in the cosmological context, is 
due to Kibble |4],[|]. Along with a precise formulation of the mechanism by which topological 
defects are produced in phase transitions Kibble made a rough estimate, using a thermo- 
dynamic equilibrium picture, of typical domain sizes. This was done by determining how 
big an ordered domain had to be in order to remain stable under thermal fluctuations to 
the disordered phase. But domain formation is inherently a non-equilibrium process and, 
therefore, one must be careful about such naive arguments. A more refined viewpoint, based 
partly on non-equilibrium ideas, has been proposed by Zurek ||. In the latter picture one 
focuses on the competition between the applied quench rate and the intrinsic relaxation rate 
of the order parameter. Defects will freeze out when the order parameter cannot adjust 
rapidly enough to follow the quench. The correlation length at this freeze-out time deter- 
mines the characteristic size of ordered domains and hence the typical separation between 
the resultant topological defects. The non-equilibrium ingredients of this picture imply that 
the domain size may differ differ significantly from the simple predictions following from 
thermal stability. Zurek and co-workers have also given numerical evidence for this picture 
using computer simulations M. 

The processes of domain formation and growth in quantum field theory at finite temper- 
atures have been investigated by Boyanovsky and collaborators using the real-time for- 



malism pQjllf . This formalism is particularly suitable for time-dependent processes like do- 
main formation for which the more popular (imaginary time) Matsubara frequency method 
is rather cumbersome to use. Applying these techniques for scalar field theories with quartic 
interactions Boyanovsky et al. |J managed to reproduce the well-known Cahn- Allen growth 
laws in the classical theory of phase-ordering |IJ. The phase transition in such systems was 
induced by adopting a sudden quench of the heat bath. Though the quench is performed 
on the heat-bath rather than on the system, Boyanovsky et al. assumed that the resulting 
effect was mocked up by taking the mass-function for the scalar fields to be a step-function 
with a change of sign. Consequently, unstable long-wavelength modes appear in the theory 
which grow exponentially with time and lead to the appearance of domains. 

Scalar field theories admitting global 0(d) symmetry admit topological defects in vari- 
ous dimensional spaces. These global topological defects contain singularities of the order 
parameter which are the zeros of the 0(d) symmetric scalar field. The density of such de- 
fects is then related to the density of the associated zeros and these can be determined, in 
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a Gaussian approximation, by the Halperin-Mazenko-Liu p2| , p!3[| formula, namely 

d/2 



n(t) 



W(0,t) 



(1.1) 



where W(x,t) is given by the diagonal equal-time two point-function ($ a (x, t)<J>fe(0, t)) = 
$abW(x, t) and the double primes denote the derivative with respect to the radial variable 
r = |r|. Once the two-point function is determined one can use (|1 . 1|) to find the defect 



density ||14|| . Though one might question the validity of the Gaussian approximation in an 
interacting system, for 0(n) symmetric field theories it can be viewed as a first order result 



in a more general scheme (15| . 
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To be precise, the use of a time-dependent mass-function is flawed as the temperature 
dependence of the masses is obtained using an equilibrium picture which is definitely not 
valid with a sudden quench. It is thus interesting to adopt a quench "slow" enough that the 
approximation to equilibrium remains valid. Though one can expect the late time behavior 
of the system to be independent of the details of the quench process, this will not be true of 
the early time behavior of the system. Our goal is to follow the system subsequent to the 
change in the mass function which we assume to be linear. Such a linear quench has recently 
been investigated in fTEfl , where it is also found that the defect densities are affected by the 
existence of a finite time quench - indeed being smaller. Though this result seems correct 
on physical grounds - its validity can be questioned due to the approximations used. On the 
other hand, the zeroth order equation of motion can be solved exactly and one can find the 
two-point function for the unstable modes, following Boyanovsky et. al. ||. We find that 
the domain sizes are larger compared to those appearing in a sudden quench, as found by 
|T5| but only during the early part of the quench and towards the end of the quench and 
thereafter the domains grow more slower in the finite-time quench. This result seems highly 
counterintuitive as a very slow quench will barely produce defects. On the other hand, the 
results of [16j show a larger domain size with increasing quench time and hence a lower 
density of defects. 

We argue that this puzzle originates from an inappropriate identification of the Goldstone 
modes which are associated with the defects - they are not affected by the quench directly. 
Indeed, we show that by adopting a polar parametrization of the global 0(2) model discussed 
in |16| - the emergent picture is consistent with the scenario proposed by Zurek. 

The outline of the paper is as follows. In section II we briefly review the results of 
Boyanovsky et. al |J, for the case of a sudden quench, to highlight the various approximations 
used. In section III we generalize to the case of a finite-time quench in their scenario showing 
that to zeroth order the domain sizes are smaller. In section IV we suggest that more physical 
results are obtained by studying the correlation functions for the polar fields. 



II. PREDICTIONS FOR A SUDDEN QUENCH 

Let us start from the well-known S0(2) symmetric 4 theory in 2+1 dimensions, which 
admits global vortices since vr 1 (S' 1 ) = Z. Most of the results below can be found in the 
paper of Boyanovsky et. al. 0. The Lagrangian for the system is given by 

L = J d 3 x[±d,cf> ■ - \m 2 (t)4> ■ <f> - -\{4> • 0) 2 ] (2.1) 
where (f> = [ } ) is a two component vector under 0(2) and we have a time- dependent 

\<P2 ) 

mass function m(t), given by 

m 2 {t) =m 2 i Q(-t) -m)Q(t). (2.2) 

Due to the 0(2) symmetry, the two-point function for this model is diagonal (0 a 0&) — 
^afc(00), where (00) is the two-point function for a single scalar field. Thus, the results for 
the 0(2) will coincide with those for a single scalar field, treated by Boyanovsky et. al. M. 
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The zeroth order equation of motion yields 

[^ + k 2 + m 2 (t)]W fc (t) = 0, (2.3) 

where U k {t) is the single particle wavefunction, as defined in the appendix. 
With the mass-function defined by fl2.2|), at times t < , Eqn. ( |2.3j ) reads, 

[^ + k 2 + m 2 ]W fe (t)=0, (2.4) 

so that one has as a solution of the form 

U k {t) = e~^\ (2.5) 

with uj 2 (k) = k 2 + m 2 , which we can treat as an "initial condition". 

In the sudden quench scenario, m 2 {t) changes instantly and at times t > 0, Eqn. ( |2 .3|) 
reads 

[^ + k 2 -mJ]W fc (t) = 0. (2.6) 

Clearly, the modes with k 2 < mj are unstable and they are responsible for the formation of 
domains fT7|| . For these modes, the solutions to Eqn.([H)|) are given by 

U k (t) = A k e w ^ + B k e~ w ^\ (2.7) 

with 



W(k) = ^m) - k 2 . (2.8) 

As we are dealing with a second order differential equation, the solutions ( |2.5| ) and 
and their first derivatives are required to match at t = , leading to 



2 y W(k) J 

B k = {A k )\ (2.9) 

where * denotes complex conjugation. For brevity, hereafter we will suppress the functional 
dependence of to and W on k. 

Given U k (t), one readily finds for t > , using 



G < r - " = \ f (2$b coth ^ l( 1 + ( w?) C08h 2Wt + i 1 - <f> 2 )] e ' k "- (2 ' 10) 

Note that in the above integral only the unstable modes are taken into account, following 

The above integral, however, includes the contributions that are already present before 
the quench i.e. for t < and hence is not appropriate for studying the growth of domains. 
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One therefore must subtract the contribution that is present at the beginning of the quench. 
For the unstable mode with wavenumber k , the growth is given by the function 



S k (t) = \U k (t)\ 2 -\U k (0)\ 2 - (2.H) 
In coordinate space this translates into 

s ^=Jwk cotH ^ )St{t)e>h " (2 - 12) 

which is obtained by replacing \U k (t)\ 2 with S k (t) in Eqn. (|A5|) . 

In the limit of high initial temperatures ( (3 — > 0) and late times, G(r, t) can be evaluated 
using a saddle point approximation, since the function ke 2Wt has a very sharp maximum at 



k max ~ yjC^) (2.14) 

This saddle point arises from the competition between the phase space factor k and the 
growth factor W(k)t. Note that this saddle point is absent in one spatial dimension. 

An expansion of the integrand around this saddle point allows us to perform the inte- 
gration at late times, leading to the Gaussian fall-off for the two point function, 

~. , m f (l + L- 2 ) . x ,e 2m f s * 2 

G(x,s) ~ L M-i=) e-w, 2.15 

4V27T V2s s 

where L = x = m^r and s = rrift. 

The domain sizes at late times can then be read off by looking at the ratio 

D(*, 5 ) = £^Wl§, (2.16) 
G(p,s) 

where we have used the fact that for s — ► oo, J (-^=) — > 1. Therefore, the domain size at a 
■ is ffiven bv 



n 

time t ^> — is given by 



ffl = J— ■ (2-17) 

V m f 

This scaling relation is the same as predicted by the Cahn-Allen equation fl|] in the 
classical theory of domain formation. The result ( |2.15| ) requires the existence of a saddle- 
point, which is absent in one dimension. It is worthwhile noting that the Cahn-Allen equation 
also fails in one dimension. In a sense, then, the saddle-point approximation is equivalent 
to the classical result. 

Now as the <3> 4 interaction and the associated back reaction effects are ignored, one has 
to be careful with this zeroth order result. In fact, the zeroth order result has the unphysical 
feature of indefinite domain growth. This shortcoming, however, cannot be overcome in a 
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perturbative framework, as emphasized by Boyanovsky et. al. || and others |18|]. One must 
treat the problem non-perturbatively using, for example, the Hartree-Fock approximation. 

To study the non-perturbative dynamics of domain formation, Boyanovsky et. al. || 
employed a time- dependent Hartree-Fock approximation. Strictly speaking such approxi- 
mations can only be justified in an ■h expansion for O(N) models. This is definitely ques- 
tionable for our case with N = 2. Despite this shortcoming, let us proceed by performing 
the decomposition 

(f>(x,t) = ^(t)+ X {x,t) (2.18) 

in the Lagrangian and retaining only the terms quadratic in the "fluctuation field" 
Note that <&(£) gives us the vacuum expectation value of the field 4>(t) at time t, i.e. 

<<p(x,t)>=$(t), <x(x,t)>=0. (2.19) 

In the Hartree-Fock approximation one assumes a condensate for the two-point function so 
that one has 

< (j) 2 (x,t) >= W(t) (2.20) 
The equation of motion for the fluctuation field is now 

[^ + k 2 + m 2 (t) + ^(t)]W fc (t)=0 . (2.21) 

We also have the condition 

$ 2 (t) = J d 2 k\U k (t)\ 2 coth(-^w). (2.22) 

The set of equations Q2.21|) and (|2.22|) have to be solved self-consistently. 

Nevertheless, one can study the qualitative behavior of the solution at late times. Note 
that for a sudden quench, the equation of motion reads, for t > 0, 

+ k 2 - m) + ^(t)]U k (t) (2.23) 

where we have chosen $ 2 (0) = 0. Since $ 2 grows with time as a result of fluctuations, the 
range of the unstable modes is reduced. At the spinodal time t = t s , defined by 



mj = -$ 2 (t s ) (2.24) 



2 

there are thus no unstable modes in the system and domain growth halts 



III. INCLUSION OF A FINITE QUENCH 

Realistically any temperature quench occurs within a finite time. How does this affect 
the above mentioned results? To answer this question, we will assume hereafter that the 
quench begins at time t = and ends at time 
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FIGURES 



m 2 (t) 




> t 



FIG. 1. m 2 (t) versus time t for a linear quench 



t = tq and it occurs linearly, so that the mass function is of the form 

m? for t < 

m 2 {t) 



2 _ t( «^ L) for < t < TQ 



mj for t > tq 
Let us now define the following quantities: 



x 



mj + m 2 



y = x 3 (jj 
z = x^W, 



(3.1) 



(3.2) 

(3.3) 
(3.4) 

where u and W are defined as in (|2.5| ) and ( |2.8| ). The solution to the zeroth order equation 
prior to the quench, for t < 0, is given by ( |2.5|) as before. During the quench period, 
t G [0, tq], however, the equation of motion at zeroth order reads, 



[^ + k 2 + ^-i W ) = o. 



(3.5) 



Let us now discuss the consequences of the finite-time quench. 
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A. During the quench 



During the quench process, the solution to Eqn. ( |3.5|) can be found in terms of the Airy 
functions, 



U k (t) = a k Ai(x a(t - u) 2 xf) +b k Bi (x z(t-u 2 x)^. 



(3.6) 



The coefficients a k , b k can be determined by matching the solutions fl2.5|) and ( |3.6|) and their 
first derivatives at the beginning of the quench , t = : 



flfc = - 

7T 



TT 



Bi'(-y*)+iyBi(-y 2 
Ai \-y 2 ) + iyAi{-y 2 



(3.7) 



where the prime denotes differentiation with respect to the argument. 

The domain growth function can be defined for the finite-time quench similar to ( p.ll|) . 
But, while in the sudden quench the unstable modes appear instantaneously, they will appear 
with a time lag, in the case of a finite-time quench. This can be understood by noting that 
in a linear quench the mass function becomes negative at time ( cf. Fig. 1 ) 



2 

mi 



TO = TQ 



m 2 + rrij J ' 



(3.8) 



with the quench beginning at t — 0. Thus the function that describes the growth of the 
unstable modes with wave-number k is given by 



S k (t-T ) = \U k (t)\ 2 -\U k (T )\'' 



(3.9) 



For an instantaneous quench there is no time lag, i.e. r = 0, and we get back the standard 
subtraction |§ in (PH . 

There is another difference between a sudden quench and a finite-time quench. Consider 
a time t , with tq < t < tq, and mass function m 2 (t) = —M 2 (t). The phase space for the 
unstable modes is given by k 2 G [0,M 2 (t)]. Note that this interval is smaller than that for 
the instantaneous quench, as 



t 



m 



f 



M 2 (t) = K + mj)[l ] > te [r , t q ]. 



(3.10) 



Hence the function describing the domain growth is given by 

rk 2 =M 2 (t) d 2 k 



G T0 (r,t) 



(27r) 2 2tu(k) 



coth( — - — )S k {t;T )e , 



(3.11) 



where the integration is performed over the unstable modes appearing at time t G [r , tq\. 
The domain size can be read off, as before, by looking at the ratio 



D T0 (x,s) = 



G T q (x, s) 
G TO (0,s)' 



(3.12) 
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In Figs. 2,3,4 a comparative plot of the domain functions for an instantaneous quench 
and a finite-time quench is given for r < t < tq. 
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0.0 1 1 1 1 1 1 1 1 1 1 — ^ 

0.0 2.0 4.0 6.0 8.0 10.0 

FIG. 2. Comparative plot of D{x,s) and D T0 {x,s) vs. x with L = .25 at s = .5 with 
sq = rrifTQ = 1, (3m f = 6. The thin line corresponds to the case with a finite-time quench 
while the thick line corresponds to the case with instantaneous quench. 
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FIG. 3. Same plot as in Fig. g but with s=.7 
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FIG. 4. Same plot as in Fig. | but with s=.9 

Note that during the early part of the quench, the domain sizes in the fmite-quench are 
larger while later the domain sizes tend to be smaller. This can be explained as a competition 
between the time lag and the reduced phase space for the unstable modes. Due to the time 
lag the domains have a shorter time period to grow while due to the smaller phase space 
the domain function has a larger spread. The latter can be understood by noting that the 
domain function is defined in position space as a Fourier transform. At early times the phase 
space is quite small and therefore washes out the effect of the initial time-lag. Domain sizes 
consequently appear larger than for a sudden quench. On the other hand, at late times, the 
available phase space for the unstable modes is comparable to that for the sudden quench. 
At late times the effect of the initial time-lag then becomes clearly visible. 

B. After the quench 

At the end of the quench, the equation of motion is given by (|2.6|) , with solutions of the 
form (fTF|). 

The coefficients A k ,B k can be determined, again, by matching the solutions (|3.6| ) and 
( f2.7| ) and their first derivatives at the point t = tq. This leads to 



Ai 



Bk 



+ h Bi(z 2 



Bi '(z 2 \ 



a k ^Ai(z 2 ) - - z Ai \z 2 ) 
a k {a%(z 2 ) + - z Ai '(z 2 ) j + b k \bi(z 2 ) + -Bi 



-tqW 



(3.13) 



The solution after the quench is thus given by 



U k {t) = a k [ Ai(z 2 ) cosh W(t - r Q ) + -Ai \z l ) sinh W(t - t q ) 
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+ b k [Bi(z 2 ) cosh W(t-r Q ) + -Bi '(z 2 ) sinh W(t - r Q )}, 



(3.14) 



where a k and b k are given by Q3.7|) . 




o.o 



0.0 2.0 4.0 6.0 8.0 10.0 



FIG. 5. Same plot as in Fig but with s = 4 



One can again use the domain growth function in (|2.11| ) but with ( |3.14|) . As the quench 
process is finished, however, all the unstable modes with < k < Mf are present after the 
quench. The domain size is given as before by ( |3.8|) . A comparative plot is given in Fig. [5|. 
This shows the domain sizes to be smaller compared with the sudden quench results, as a 
consequence of their slower growth during the quench period. 

This result, however, seems to be in contradiction with the physical picture in which a 
very slow quench should not lead to appreciable domain formation. This puzzle points out 
that we have to take into account the effects of the back-reaction. 



C. The role of Back- Reaction 

Let us now incorporate the back-reaction into the picture with a finite-time quench. The 
equation of motion in the Hartree-Fock theory, during the quench period, would read as 

[^ + k 2 + m, 2 + ^(t)-^(t). (3.15) 

Solving this equation is difficult. To gain some qualitative understanding, let us again 
consider an adiabatic quench. Note that as $(t = 0) = there is not much change in 
the initial behavior of system compared to the zeroth order case. As $(t) grows, however, 
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one has $ 2 (£) > and the "back- reaction" competes with the decreasing mass-function, 
producing a slower quench. Recall that domain formation is associated with the appearance 
of unstable modes and due to the slower quench they will appear even later compared to 
the zeroth order case. This is definitely true for small values of A, so that the mass function 
can pass through zero at late times. By our earlier arguments, one can see that the domains 
formed subsequently will be larger again. On the other hand, for larger values of A, if 
$ grows fast, the effective mass-function might not become negative at all during a weak 
quench and there will be no domain formation. Thus, if defects are to form one would need 
a quench rate faster than the growth rate for the field $, as envisaged in the picture of Zurek 



IV. THE QUENCH ON THE GOLDSTONE MODES 

The Halperin-Mazenko-Liu defect density formula involves the complete two-point func- 
tion whereas we have been dealing with a restricted two-point function which only involves 
the long-wavelength unstable modes. On the other hand, defects are localized objects and 
thus they carry short-wavelength modes with them. In fact, the 0(2) global defect is asso- 
ciated with a phase variable which has no apparent cutoff associated with it, as we shall see 
below. 

In our above treatment of the 0(2) symmetric scalar model, we have used a Cartesian 
parametrization for the fields. The zeroth order equations of motion for the two fields 
accordingly decouple and appear as linear equations. The phase degree of freedom (the 
Goldstone mode of the theory) nevertheless interacts with the radial mode. The defect 
configurations in the 0(2) model are associated with windings of this Goldstone phase. 
It can therefore be misleading to use the Cartesian parametrization when studying defect 
formation. Let us instead rewrite the 0(2) model in terms of a complex field $ (polar 
decomposition): 

$ = 1 + 20 2 = Fe ie . (4.1) 
In terms of F and O, the Lagrangian now reads 

L=- [d^Fd^F + F 2 d tl Qd^Q] - ^m 2 (t)F 2 - ^XF\ (4.2) 

exhibiting the Goldstone nature of the phase O and its coupling to the radial mode F. 
The equations of motion following from the Lagrangian (|4.2|) are given by: 



d^F + [m 2 {t) + AF 2 - (0 M ed"e)] F = (4.3) 

dJF 2 d»Q) = (4.4) 



Note that in (O) the term involving the Goldstone modes appears with the opposite sign 



to the time-dependent mass function m 2 (t), in contrast to the case with back-reaction. Note 
also that the quench does not affect the Goldstone modes directly. 

We will attempt to follow the dynamics of this system by using a Born-Oppenheimer 
like approximation, treating O as a fast variable and F as a slow variable. 
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If we ignore the back- reaction due to the field G, the VEV for the radial field F is given 
at any instant during the quench by 



AF 2 m-i° {t = 0) UK\ 

A K } ~ \ -m 2 (t) + (fyea^e) (t > o) 1 0> 

Using (|4.5|) in the second equation of (|4.4|) one gets — 

A<9 M (F 2 ^6) 

= <9 M [(m 2 (t) + (d u ed u e))d"] 6 = (4.6) 

One can attempt to simplify the dynamics by ignoring the (d u Od u O) term in the above 
equation and solve instead the equation, 

d,(m 2 (t)d»e) =0, (4.7) 

where we have dropped the extra term from the equation of motion as it is less relevant at 
long-wavelength. Equation (|4.7|) can be expanded as 



e + ^-(^m 2 )e-v 2 e = o. (4.8) 

m 2 at 

The second term can be interpreted as a "friction" term which changes sign from positive 
to negative as the mass function changes sign. Accordingly, there will be production of 
Goldstone particles which will feed into the dynamics of the F field, as in ( |4.3|) . This is 
similar to the picture in ||19|| . Since there is no intrinsic mass scale for 9 field one needs to 
consider all frequencies including the high ones. This is also necessary as the defects are 
localized objects which can only be probed with high frequencies. 

Note that when m 2 (t) — > , only the second term which involves the first derivative in 
time dominates, 

2m6 ^ for m -> 0. (4.9) 

In the linear quench this shows that the Goldstone modes become time-independent imme- 
diately prior to the appearance of the unstable modes. Thus the defects get disentangled 
from the evolution of the field F, as in the picture proposed by Zurek ||. 
An adiabatic solution to Eqn. (|4.8|) is given by 



where 



0(x,t) = J d 2 kduj[6(uj,k) + c.c], 

6(u,k) = e lk -*e~ iuJt (4.10) 



^ 4m 4 i 2 ^ 2m?x ^ ' 



in the adiabatic approximation, with x = m i+ m i ■ For the static situation, one has x = oo 

and one gets back the standard free-particle. On the other hand, when m 2 — > 0, we find 
uj — > 0, — zoo leading to the static situation discussed above. 



13 



Though qualitatively the above scenario supports the picture of defect formation, there 
is definitely much work needed to get quantitative predictions. Specifically, the domain sizes 
can be found from the two-point function for the F field. Because of the coupling to the 
Goldstone modes the evaluation of this two-point function is quite involved. One alternative 
would be to integrate out the Goldstone field ( after a Gaussian integration) leading to 
an effective action involving the F field only 

L e// (F) = L -^Tr ln(^F 2 ^) 

= L - -Tr ln(F 2 a 2 ) - -Tr ln(l + + . . . . (4.12) 

Note that in the second line we have expanded the logarithm. Note that the operators 
F 2 d 2 and d tl F 2 d >1 do not commute and therefore in the above expansion (using the Baker- 
Campbell-Hausdorff formula) there are other terms which are represented by the ellipsis. 
The second term basically cancels the dependence on F in the functional measure in the 
path integral. Though naively it seems that the third term in ( 4.12 ) suffers from an infrared 



divergence, the finite temperature of the heat bath provides us with an infrared cutoff. In 
light of this one can attempt to expand the logarithm in the third term to get an effective 
action for F which takes into account of the "backreaction" of the Goldstone modes. In 
practice this is rather complicated as the the functional determinant is time-dependent (due 
to the time-dependence of the cutoff). For an initial high temperature, the logarithm can 
be expanded as ~ is small. During the later stages of the quench, however, 

i 

this may no longer be true. It might be more appropriate, therefore, to use the recently 
developed formalism of the non-equilibrium effective action []2{J since it generates time- 
dependent correlation functions. 



V. CONCLUSIONS AND REMARKS 

In this letter, we have examined the formation of defects in phase transitions induced 
by a finite-time quench on 0(2) symmetric scalar field theories. We find that adopting the 
treatment of Boyanovsky et. al. naively with a Cartesian parametrization leads, to the 
zeroth order, to smaller defect densities during the early quench period while higher defect 
densities for the later phase. This shows that the problem of the back-reaction needs to be 
studied properly. 

We also propose that a more physical picture will be found by adopting the polar 
parametrization of the group manifold. This makes explicit the interaction of the true 
Goldstone phase field with the time-dependent radial field. In the polar parametrization the 
theory becomes interacting even in the absence of the $ 4 interaction. This makes the anal- 
ysis of the model much harder. For a finite time quench, however, such a parametrization 
suggests qualitative support for the Zurek picture of defect formation. This will be valid for 
the 0(n) symmetric scalar models as well. One can attempt to solve the time-dependent 
model numerically. We hope to report on this in the near future. 
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APPENDIX: REAL-TIME FINITE TEMPERATURE FIELD THEORY 

The real time formalism ( also known as the Schwinger-Keldysh formalism) allows one to 
compute in-in matrix elements, and hence is useful for computing time-dependent averages. 
For a review of finite temperature field theory in this formalism see, for example, |TD[. Here 
we only quote the formulae needed for our purposes. 

In this paper, we will be interested in the equal-time two-point Green function 

G{x, t) = < (f>(x, t)(f){0, t)> = Tr [p <p{x, t)<p{0, t)}. (Al) 
Here po is the thermal density matrix describing the initial condition, 

e -0Hi 

with Hi the Hamiltonian describing the system initially. 

For a bosonic field in D spatial dimensions with a mass m we have the standard expan- 
sion: 

where u(k) = (k 2 + m 2 )^ and Uk(t) is the one-particle wavefunction in momentum space. 
The particle creation and destruction operators satisfy 

[a k , al] = {2ir) D 2u{k)5 D {k - k'). (A4) 

Substituting (|A3|) in ( [AID one gets, 
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